Analysis of Exactness of the Mean-Field Theory in Long-Range Interacting Systems 
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Relationships between general long-range interacting classical systems on a lattice and the cor- 
responding mean-field models (infinitely long-range interacting models) are investigated. We study 
systems in arbitrary dimension d for periodic boundary conditions and focus on the free energy for 
fixed value of the total magnetization. As a result, it is shown that the equilibrium free energy of the 
long-range interacting systems are exactly the same as that of the corresponding mean-field models 
(exactness of the mean-field theory). Moreover, the mean-field metastable states can be also pre- 
served in general long-range interacting systems. It is found that in the case that the magnetization 
is conserved, the mean-field theory does not give correct property in some parameter region. 



Long-range interactions appear in several fields [1] : as- 
trophysics, plasmas, 2D hydrodynamics, and so on. Re- 
cently it was reported that a model of spin-crossover 
materials has an effective long-range interaction among 
molecules due to the coupling to the local lattice distor- 
tion [11]. It has been shown that statistical mechanics of 
long-range interacting systems exhibits several peculiar 
features: negative specific heat [2], long-lived metastable 
states [3], and ensemble inequivalence [4]. As simpli- 
fied models to study these unfamiliar features, the so 
called mean-field models (infinite-range interacting mod- 
els) have been often adopted. From the statistical me- 
chanical point of view, one expects that at least some 
qualitative natures of systems with slow decaying long- 
range interactions are captured by the analysis of corre- 
sponding mean-field models [1]. Indeed, some evidences 
for the fact that the mean-field model gives an exact de- 
scription of equilibrium states (exactness of the mean- 
field theory) in the models with l/r Q -type long-range in- 
teraction (0 < a < d, where d is the spatial dimension- 
ality) have been obtained by the studies of specific mod- 
els [5-8]. There is a conjecture that the equilibrium prop- 
erties of long-range interacting spin systems are exactly 
the same as those of the corresponding mean-field mod- 
els [5]. In this letter, we address the question whether 
exactness of the mean-field theory holds in general. 

We consider a lattice system. Each lattice site i has 
a "spin " variable <r,-. In the Ising model, a.- L = ±1, but 
we do not restrict the variables {a{\ to the Ising spins. 
They can also take continuous or vector variables. For 
example, cr, = (cos 0, sin 0), < < 2ir corresponds to 
the XY model and Gi = e a {a = I, . . . ,q) corresponds 
to the g-state Potts model, where {e a } are unit vectors 
satisfying e a ■ ep = 8 a ^. We investigate the system of 
the following type Hamiltonian: 
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case J > 0. The uniform magnetic field is denoted by H . 
We define the distance |r,-— Tj\ between the lattice points 
i and j as the shortest distance of these lattice points in 
periodic boundary conditions. We consider the follow- 
ing two types of long-range interactions: the power-law 
interaction K(r) oc l/r a , < a < d, and the Kac inter- 
action [9] K(r) oc 7 d 0(7r). Here, (f>(x) is assumed to be 
non-negative 4>{x) > and integrable J d d x<fi(x) < +oo. 
Moreover, we assume that there is a positive and de- 
creasing function tfj(x) such that |<^'(x)| < ij){x) and 
J d d xtp(x) < +oo. These assumptions are necessary to 
justify the coarse-graining of the Hamiltonian (to be ex- 
plained below). A typical example of the Kac interaction 
is the exponential form, K(r) oc r ) d e~ ir . In this case, 
4>{x) = ip(x) = e~ x . 

We take the limit 7 — > in the Kac interaction. We 
consider the two limiting procedures: the van der Waals 
limit, 7—^0 after L — > 00 [10] and the long-range limit, 
7 — > with 7L = const. The former limit corresponds 
to the situation that the interaction range 7 -1 is much 
longer than the microscopic length scale (the lattice in- 
terval) but much shorter than the system size L. The 
latter limit corresponds to the situation that the inter- 
action range is comparable with the system size. These 
two limits lead to different kind of behavior. 

Note that the thermodynamic limit does not exist for 
long-range interactions in the usual sense. To restore the 
thermodynamic limit, we adopt the Kac prescription [9], 
,q K (|rj|) = 1, where the interaction is normalized 
and depends on the system size. 

We focus on the free energy restricted to a fixed mag- 
netization m, 



exp(-0F(m,T,H))= £\ 
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Here, K(\ri — rj\) is the interaction potential and J is 
the interaction strength. We consider the ferromagnetic 



which corresponds to the Landau free energy, and we call 
F(m,T, H) merely "free energy" hereafter. The param- 
eter /3 = 1/T is the inverse temperature and we set the 
Boltzmann constant to unity. The symbol 6(a, b) is the 
Kronecker delta. For conserved systems (i.e. systems 
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whose magnetization is fixed), m is a fixed parameter 
and the term —Hm is just a constant. Then, F(m, T, H) 
gives the equilibrium free energy. 

In non- conserved systems (i.e. systems whose mag- 
netization is not fixed), the equilibrium free energy is 
given by the minimum value of F(m, T, H) with respect 
to m. Local but not global minima are interpreted as 
metastable states. In preceding works, exactness of the 
mean-field theory is investigated only for non-conserved 
systems [5-8]. The merit to consider F(m, T, H) is that 
we can treat both conserved and non-conserved systems. 
It should be noted that the ensemble with fixed magneti- 
zation and that with fixed magnetic field are essentially 
different in long-range interacting systems. Moreover, it 
allows us to discuss metastable states as local minimum 
points of the free energy. 

In the long-range interacting systems, one expects that 
only long wavelength modes play important roles for 
macroscopic behavior. In fact, it is possible to perform 
coarse-graining exactly for long-range interacting mod- 
els [6]. Let us divide the lattice system into blocks of 
the linear dimension I. The number of blocks is (L/l) d , 
where L denotes the system size, and each block has l d 
sites. We introduce a local coarse-grained variable m k as 
m k = jj J2ieB k °i in each block B k,k = l,2,..., (L/l) d . 
We define the position x k = r^/L, where is the cen- 
tral position of a block B]~. We also define m(xk) = mk- 
We take the limit L — > oo, I — > oo with l/L — > 
(continuous limit). In this limit, x k becomes a con- 
tinuous variable x. For long-range interacting models, 
in a proper choice of l(L) as a function of the system 
size L such that lim£_>. 00 l(L)/L — 0, the Hamiltonian 
is written only by m(x) in the thermodynamic limit [6]: 
H = H[m(x)] + o(N), where 

T-L[m{x)] = — / d d x I d d yU{x - y)m{x)m{y) 

-NH [ d d xm(x). (3) 
Jo 

Here the interaction potential U (x) is given by 

U(x) = lim L d K(Lx). (4) 

L— 7-oc 

The integration in (3) means d d x = dx\ . . . dxd- 
The fact that the Hamiltonian is written only by the 
coarse-grained magnetization m(x) means that for a 
fixed system size L (it is very large) , there is some length 
scale 1{L) in which the magnetization is uniform, though 
the whole system may be inhomogcneous. For example, 
in the case of the power-law interaction K(r) ~ l/r a , 
this length scale is given by l(L) ~ L x ~ a l d . While the 
size of a block l(L) is itself macroscopic when L — > oo, the 
blocks are not independent each other. It means that we 
cannot divide a macroscopic system into two macroscopic 
subsystems without any macroscopic change. 



Performing the Fourier expansion in Eq. (3), we obtain 
the following expression: 

n 

where m(x) = ^ n m n e 2lTm ' x and U n = 
f Q d d xll(x) cos(27rn • x) with n e Z d . Here, the 
Kac prescription implies that U = J Q d d xU(x) = 1. 
We call {U n } the interaction eigenvalues. Dividing the 
Hamiltonian into two parts: the mode of n = and the 
others, we have 

H=(-^mt-NHm )-^fj2 U ^\ 2 - ( 6 ) 

We call the first term of the above equation the mean- 
field model corresponding to (1). 

The form of the coarse-grained potential U(x) is given 
by U(x) = A/x a for the power-law interaction, U(x) = 
S(x) for the Kac interaction with the van der Waals limit, 
and U (x) = B-fQ<f>(jox) for the Kac interaction with the 
long-range limit. The constants A and B are determined 
by the normalization condition Uo = 1. Note that U n = 1 
for any n for the Kac interaction with the van der Waals 
limit. It is known that the free energy is given by the 
mean-field theory with the Maxwell construction in the 
van der Waals limit [10]. Therefore, the free energy does 
not depend on its detailed interaction form in the van der 
Waals limit. On the other hand, for the long-range limit 
or the power-law interaction, whose interaction range is 
comparable with the system size, the free energy docs 
depend on the interaction form, as we see below. The 
modes with < U n < 1 play an important role. 

Now, we evaluate the free energy F(m,T, H). Since 
the coarse-graining can be exactly performed, the Hamil- 
tonian H can be replaced by H[m(x)]. The summa- 
tion of (2) is divided into two parts, f„ Vm(x) 

and E^} with fixed m (xY namely, the profile of m(x) 
and the configurations inside the blocks. Note that 
the Hamiltonian depends only on m(x). Writing the 
number of states with the fixed magnetization m by 
W = exp(S(m)), where S(m) is the entropy, we have 

E {CTl} with fixed m(x) 1 = c *p(Jo S{m(x))d d x y There- 
fore, we obtain 

e -0F{m,T,H) = I V m (x)e-^l?n( X ).T.H] i (7 ) 

J iriQ—m 

where J 7 [m(x),T,H] = H[m(x)} - T J* S(m(x))d d x is 
called the free energy functional. By using the saddle- 
point method, the free energy is given by the following 
minimization problem: 

F(m,T,H)= mm T[m{x),T,H}. (8) 

{m(x) |mo— m\ 
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The free energy functional is the same as the mean- 
field free energy Tmf when m{x) = m, therefore, the 
upper-bound of the free energy is obtained, F(m, T, H) < 
F MF {m,T,H). We also found that F(m, T, H) < 
F MF {m,T,H) when d 2 F MF {m,T/U max , H)/dm 2 < 0, 
that is, F(m, T, H) does not agree with Fmf(w, T, H). It 
is because the uniform magnetization profile m(x) = m 
gives the local maximum value of the free energy func- 
tional when <9 2 Fmf(«t-, T/[/ max , H)/dm 2 < 0. Here, 
U m ax < 1 is the largest interaction eigenvalue except for 
U a = 1. 

Taking into account of the relations X^n^o Un\m n \ 2 ^ 

PmaxEn/O l™™| 2 and Jo m ( x ) 2 d d X = m 2 +J2n^O \™n\ 2 , 

the free energy functional is evaluated as 



(a) 

F MF (m,T/[/ max ,0) 




"T MF (m,T,0) 




FIG. 1: (Color online) An illustrative explanation of the re- 
gions A, B, and C in the Ising model at H = 0. The blue lines 
represent the mean-field free energy at the temperature T. 
The red lines represent the mean-field free energy at the effec- 
tive temperature T/Umax- (a) The case of T/t/ ma x > Tc > T. 
(b) The case of T c > T/U max > T. 



T[m(x),T,H] >F MF (m,T,H) - U n 



F M F{m,T efi ,H) 



d d xF MF (m(x),T eS ,H) 



(9) 



where T c g = T/U max . Here the following equality holds: 



mm 

{m(x)\rho— m} 



d d xF MF {m(x),T eS ,H) 



= CE{F MF (m,T cS ,H)}, (10) 



where CE means the convex envelope. The convex enve- 
lope of a function g(x) is defined as the maximum convex 
function not exceeding g{x). Thus, it is concluded that 
F(m, T, H) > F M F(m,T,H) -U max AF M F(m,T eB ,H), 
where we define AFmf = -Fmf — CE{Tmf}. It gives 
the lower bound of the free energy. Finally, we obtain 

F MF (m, T, H) - U max AF MF (m, T/U max , H) 
< F{m, T, H) < F MF (m, T, H). (II) 

Only the inequality U n < C/ max is used to derive the 
lower bound. In the Kac interaction with the van der 
Waals limit, U n = 1 for any n, and hence the lower 
bound in Eq. (11) is realized and thus, F{m,T,H) = 
CE{Fmf(™, T, H)} holds. This result is consistent 
with [10]. 

According to Eq. (11), the parameter region (m, T, H) 
is classified into the following three regions: 

Region A: the region where AFmf (jn, ^ T , H^j = 0. 
In this region, the mean-field model gives the exact 
free energy, F(m, T, H) = F MF (m,T,H). 

Region B: the region where AFmf (m, f7 T , H^j > 

and £t^mf (m, x£^,h\ > 0. In this region, 
the uniform configuration is locally stable but 
not necesarily globally stable and F(m, T, H) < 
F MF (m,T, H). 



Region C: the region where -^-sFmp (rn, ^ , H^j < 
0. In this region, the uniform configuration is un- 
stable and the free energy is not given by the mean- 
field free energy, F(m, T, H) < F U F(m,T,H). 

Notice that this classification is determined only by 
Fmf and C/ max . We can calculate both Fmf and U ma , x for 
specific models. Hence, we can specify these three regions 
explicitly. In Fig. 1, a typical example of the regions A, 
B, and C are depicted. In the case T e g > Tq > T where 
Tc is the mean- field critical temperature (Fig. 1 (a)), 
Fmf(w, T e ff, H) is a convex function, and thus the whole 
m is in the region A. On the other hand, in the case Tc > 
Toff > T (Fig. 1(b)), the region where F M F(m, T cS , H) is 
concave exists, and it is the region C. The region between 
A and C is the region B. In the region A and a part of B, 
F(m, T, H) = F MF {m,T,H) holds. We call this region 
"the MF region" . On the other hand, in the region C and 
the other part of B, F(m, T, H) ^ F MF {m,T,H). This 
region is "the non-MF region" . 

In conserved systems, F cq (m,T) = F(m,T,H = 0) 
and the equilibrium property of the long-range interact- 
ing system is exactly the same as that of the correspond- 
ing mean-field model in the MF region, but they are not 
the same in the non-MF region. In the MF region, the 
magnetization is uniform (see Fig. 2 (a)). In the non-MF 
region, inhomogeneity appears, which is demonstrated 
in Fig. 2 (b), because the deviation from the mean- field 
model is due to the non-zero wavenumber components 
rhn^o of the magnetization. In fact, a kind of phase tran- 
sition occurs in the boundary of the MF and the non-MF 
region, and the cluster is formed at equilibrium in the 
non-MF region. However, this clustering should not be 
understood by the phase separation. In long-range inter- 
acting systems, a part of the system is not independent of 
the other part of the system, and we cannot separate the 
system into two parts without any macroscopic change. 
Anyway, in conserved systems, the long-range interact- 
ing model is not fully describable by the mean-field model 
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(a) 



(b) 




FIG. 2: Typical equilibrium configurations in the 2D lattice 
gas model (conserved Ising model) with 1/r type long-range 
interaction. The black points correspond to the occupied state 
Oi = +1. The parameters are set to be m = 0.5, J = 2 and 
L — 100. (a) The case of T = 0.6. These parameters belong 
to the MF region, (b) The case of T = 0.4. These parameters 
belong to the non-MF region. 



(a) 



(b) 
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FIG. 3: The comparison between periodic and free boundary 
conditions by the Monte Carlo calculations. The system size 
is L = 200. (a) Equilibrium magnetizations in non-conserved 
systems, (b) Equilibrium internal energies in conserved sys- 
tems at m = 0. 



even in equilibrium. 

In non-conserved systems, the equilibrium magne- 
tization m oq is determined by mhi m [F(m, T, H)] — 
F(m eq ,T,H). Because Ai*MF(m eq , T e s, H) = and all 
the equilibrium states belong to the region A, it is con- 
cluded that the equilibrium property is always equivalent 
to the mean-field model, which was conjectured by Can- 
nas et al. [5]. Moreover, our result implies that many 
(but not all) mean-field metastable states are also pre- 
served in general long-range interacting systems because 
local minima of the free energy are in the region A or B. 
Indeed, we have checked by numerical calculations that 
in the Glauber dynamics, relaxation times from the local 
minima of the free energy to the equilibrium states are 
identical to those of the corresponding mean-field model 
(not shown) . To study some universal dynamical aspects 
is an important future problem. 

Finally, we comment on boundary conditions. We im- 
posed the periodic ones in this letter, but boundary ef- 
fects will be important for long-range interacting sys- 
tems. Let us redefine the interaction eigenvalues for free 
boundary. In this case, the interaction eigenvalues are 
not given by the Fourier coefficients. In general, we re- 
gard Uij = K(\ri — Tj\) as a matrix and denote its eigen- 
values by {uk} (uq > u% > U2 > • • • )• We call Uk as 



an interaction eigenvalue and normalize the interaction 
to set uq = 1. Regarding ui as U max , we have found by 
Monte Carlo calculations that the above redefinition of 
interaction eigenvalues allow us to extend our result to 
the free boundary case. For example, see Fig. 3. We con- 
sider the two dimensional Ising model (er = ±1) on the 
square lattice with the 1/r type interaction. In this case, 
U m ax = 0.31 for periodic boundary and J7 max = 0.42 for 
free boundary. In Fig. 3 (a), the equilibrium magnetiza- 
tions in non-conserved systems are depicted. Our result 
states that the mean-field theory is exact and the phase 
transition occurs at T/J = 1. In Fig. 3 (b), the equi- 
librium energies in conserved systems (we set m = 0) 
are depicted. Our result states that when T/J > C/ max 
(the region A), the system is uniform and the energy is 
constant, but when T/J < J7 max (the region C), the in- 
homogeneity grows. These numerical results are consis- 
tent with our result for both boundary conditions. The 
boundary effects will be examined in more detail else- 
where. 

In conclusion, we have studied the robustness of the 
results of the mean-field model against the interaction 
forms. We revealed that the results of the mean-field 
model are fully reliable in a wide parameter region (the 
MF region), but there is the other region (the non-MF 
region) where the mean-field model cannot describe long- 
range interacting systems. Properties of states in the 
non-MF region are not obvious yet. It is a future problem 
to study them. 
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